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Abstract.  A  “Carathcodory-Fcjcr  method”  is  presented  for  near- best  real  ra¬ 
tional  approximation  on  intervals,  based  on  the  eigenvalue  (or  singular  value)  analysis 
of  a  Ilankcl  matrix  of  Chebyshev  coefficients.  In  approximation  of  a  smooth  func¬ 
tion  F  ,  the  CP  approximant  Rc *  frcquentljKdilfcrs  from  the  best  approximation 
II *  by  only  one  part  in  millions  or  billions.  To  account  for  this  we  show  here  un¬ 
der  weak  assumptions  that  if  F  is  approximated  on  [— c,t]  ,  then  as  e  — ►  0  ,  \\F  — 
/C'll  =  0(<m+n'M)  while  || Rcl  —  7?*||  =  0(e3m+2n+3)  [n  contrast,  the  latter  figure 
would  be  0(cm4’n+J)  for  the  Chebyshev  economization  approximant  of  Maehly  or  the 
Chebyshev- Pad 6  approximant  of  Gragg.  It  follows  that  as  c  — ►  0  ,  best  approxima¬ 
tion  error  curves  approach  the  real  parts  of  m  +  n  4-  1  -winding  Blaschke  products 
to  within  0(c3m+2n+3) .  Numerical  examples  are  given,  including  applications  to 
ex  on  ( —  1 ,  l]  and  e“x  on  [0,  oo) .  For  the  latter  problem  we  conjecture  that  the 
errors  in  (n,  n)  approximation  decrease  with  each  n  by  a  ratio  approaching  a  fixed 
constant  9.28903.... 
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0.  INTRODUCTION 


The  purpose  of  this  paper  is  to  describe  and  analyze  a  new  analytical  method  for  near¬ 
best  rational  Chebyshev  approximation  on  an  interval,  which  we  call  the  Caratheodory- 
Fejer  method,  that  is  based  on  an  eigenvalue  analysis  of  a  Ilankel  matrix  of  Chebyshev 
coefficients.  What  is  remarkable  about  the  CF  method  is  its  extraordinary  degree  of 
optimality  in  approximating  many  smooth  functions.  Let  R*  be  the  best  (Chebyshev) 
approximation  of  rational  type  (m,n)  on  [—1,1]  to  a  continuous  function  F{x) ,  let  Rc ? 
be  the  corresponding  CF  approximation,  and  let  Em  =;  \\F  —  /?*||oo  and  Ec *  =  ||F  — 
Rcf\\oo  be  the  associated  errors.  Then  Table  1  shows  how  close  Ec*  and  E *  turn  out  to 
be  for  the  case  F(x)  =  e*  . 


cf 

(m,n)  £*  E  -  E*  (approx.) 


(0,0) 

1.1' 32 

IQ"4 

(1,1) 

2.0970(-2) 

IQ'6 

(2,2) 

8.6900(-5) 

io~12 

(3,3) 

1.5507(-7) 

10-20 

(4,4) 

1.538K-10) 

< 10-27 

Table  1.  Errors  in  best  and  CF  approximation  of  eX  on  (jl,l] 
by  rational  functions  of  type  (n,n)  ,  0  <.  n  <  4  . 


Such  extremely  strong  agreement  demands  explanation.  In  this  paper  it  is  shown  that  if  a 
smooth  function  F(x)  satisfying  a  simple  normality  condition  is  approximated  on  [— e,c] , 
then  || Rc*  —  72* ||  =  0(c3m+2n+3)  as  e  — ►  0  (Thm.  6).  As  a  corollary  it  is  also  shown 
that  to  the  same  order  as  c  — ►  0  ,  F  —  Rm  equals  the  real  part  of  an  (m  +  n  +  1)  -winding 
Blaschkc  product  on  the  complex  unit  circle  (Thm.  7). 

Our  research  on  CF  approximation  began  with  the  study  of  error  curves  in  approxima¬ 
tion  of  analytic  functions  on  the  complex  unit  disk.  If  rm[z)  is  the  best  approximation 
of  type  (m,n)  to  f{z)  on  \z\  <  1 ,  it  turns  out  that  the  error  curve  (/  —  =  I) 
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often  approximates  extremely  closely  a  perfect  circle  about  the  origin  of  winding  number 
m  +  n  +  l  .  This  phenomenon  for  polynomial  approximation  was  discussed  by  Trefethen 
in  [40],  where  by  an  analysis  based  on  the  Caratheodory-Fejer  theorem  [2,5],  it  was  shown 
that  in  best  approximation  on  the  disk  |z|  <  e  ,  the  error  curve  is  circular  to  0(e2m+3)  as 
€  — ►  0  .  By  means  of  an  extension  of  the  CF  theorem  due  originally  to  Takagi  [19,37]  and 
generalized  by  Adamjan,  Arov,  and  Krein  [1],  this  result  was  extended  to  0(c2m+2n+3) 
for  rational  approximation  in  [41].  At  the  same  time  Gutknecht  found  that  the  CF  tech¬ 
nique  could  be  transplanted  by  the  Joukowski  map  x  =  \{z  +  z~l)  from  |z|  =  1  to  x  6 
[—1,1],  and  the  resulting  real  CF  method  was  analyzed  for  polynomial  approximation  in 
[21].  The  present  paper  completes  this  series  by  presenting  and  analyzing  asymptotically 
a  CF  method  for  real,  rational  approximation. 

The  Joukowski  transplantation  has  previously  been  applied  for  near-best  real  ap¬ 
proximation  by  Frankel,  Gragg,  and  Johnson  [14,17],  who  derived  a  Chcbyshev-Pade  ap¬ 
proximation  on  [—1,1]  based  on  Pade  expansions  at  x  =  0.  This  Chebyshev-Pade  ap¬ 
proximation  is  related  to>  but  not  the  same  as,  the  earlier  rational  economization  fraction 
of  Maehly  [6,  p.  178;  17].  Our  fraction  might  be  called  the  Chcbyshev-CF  approximant, 
for  it  fits  directly  into  the  framework  of  Gragg  and  his  colleagues.  Indeed,  correspond¬ 
ing  to  their  Fourier- Pade  and  Laurent- Pade  approximations,  one  can  develop  Fourier-CF 
and  Laurent-CF  approximations  for  real  periodic  and  complex  meromorphic  functions, 
respectively.  In  general  the  CF  approximations  will  be  more  complicated  but,  for  smooth 
functions,  much  closer  to  optimal. 

The  history  of  ideas  connected  with  the  real,  rational  CF  method  is  long  and  confused. 
Eigenvalues  of  Ilankel  matrices  were  used  half  a  century  ago  for  estimating  the  error  of 
the  best  approximation  and  for  solving  certain  special  problems  exactly  by  Bernstein, 
Achieser,  and  Mirakyan;  see  [31,  p.  106]  and  [2,  App.  D]  for  references.  The  use  of 
such  a  device  for  near-best  approximation  was  apparently  first  proposed  by  Darlington 
in  1970  [10]  for  the  real  polynomial  case,  and  the  .first  (and  only  previous)  extension  to 
rational  near-best  approximation  is  due  to  Lam  and  D.  Elliott  in  1972  3  [12,27,28].  The 
connection  between  the  CF  method  and  approximation  on  the  disk  was  first  pointed  out 
in  the  excellent  dissertation  of  Ilollenhorst  [2.5]  (for  the  polynomial  case),  and  this  was 
also  the  first  work  to  contain  error  estimates.  Further  related  contributions  have  also 
been  made  by  C.  Clenshaw,  G.  H.  Elliott  [13],  A.  Reddy,  and  A.  Talbot  [38,39],  One  of 
our  own  contributions  in  previous  papers  [21,40,41]  has  been  to  connect  CF  methods  with 
the  Caratheodory-Fejer  theory  and  the  related  results  of  Takagi  and  Adamjan,  Arov,  and 
Krein.  This  makes  it  possible  to  fill  various  theoretical  gaps.  A  second  innovation  in  our 
papers  has  been  that  by  means  of  arguments  related  to  strong  uniqueness,  we  apply  results 
of  the  CF  method  to  get  estimates  on  the  behavior  of  the  best  approximation  itself.  In 
our  view,  the  CF  idea  is  not  just  a  method  for  generating  near-best  approximations,  but  a 
theory  that  should  reveal  hitherto  unrecognized  properties  of  real  and  complex  Chebyshcv 
approximation  and  of  the  relation  between  the  two. 

The  present  paper  also  differs  from  the  previous  work  of  Lam  and  Elliott  in  two 
practical  ways.  First,  our  method  applies  for  arbitrary  m  and  n ,  rather  than  just  m  > 
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n  ;  the  idea  behind  tills  extension  is  derived  from  an  example  by  Talbot  [38].  Second,  its 
asymptotic  order  of  accuracy  on  small  intervals  is  0(c3m+2n+3)  rather  than  (it  appears) 

0^2m+2n+3  j  ^ 

Our  arguments  proceed  in  two  steps.  First,  one  shows  that  the  CF  method  yields  an 
error  curve  that  nearly  cquioscillates  (real  case)  or  is  nearly  circular  (complex  case);  this 
implies  Ec ?  «  E *  .  Second,  one  shows  by  an  argument  related  to  strong  uniqueness  or 
Lipschitz  continuity  of  best  approximations  that  this  behavior  further  implies  Itcf  ^  R*  . 
All  of  our  estimates  are  asymptotic,  pertaining  only  to  the  interval  [ — c,  c]  in  the  limit 
e  — ►  0  .  (Rquivalontly,  one  could  consider  increasingly  smooth  functions  F(fx)  on  the 
fixed  interval  [ — i,  t] .)  This  is  effectively  the  same  limit  considered  in  the  past  in  various 
papers  on  Chebyshev  approximation,  notably  [29,34]  for  real  polynomial,  [7,20,43]  for  real 
rational,  [32]  for  complex  polynomial,  and  [42]  for  complex  rational  approximation.  Where 
these  papers  obtain  one  term  of  an  asymptotic  expansion  of  the  best  approximation  (two, 
in  the  case  of  [34]),  the  CF  method  gets  many. 

The  CF  method  and  related  ideas  are  currently  attracting  much  attention  in  the 
theories  of  digital  filtering,  control,  and  linear  systems.  This  work  has  been  mainly 
stimulated  by  the  papers  of  Adamjan,  Arov,  and  Krein  [1],  and  is  being  carried  out  by 
(among  others)  M.  Bcttaycb,  A.  Bultheel,  P.  Dewilde,  Y.  Genin,  S.  Kung,  and  L.  Silverman. 
See  the  book  by  Kailath  [26]  for  some  references,  and  also  [20]. 
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1.  TIHi  CARATIII30D0RY-FI3JER  METIIOP 

Let  the  unit  disk  and  circle  be  denoted  D  =  {z  :  |z|  <  1}  and  dD  =  {z  \  |z|  =  1}  , 
let  dD+  be  the  upper  semicircle  dDf \{z  :  Imz  >  0} ,  and  let  /«  =  [— c,  c] .  Let  ||  •  ||/,  , 
II '  llao  >  etc.  be  the  corresponding  supremum  norms,  but  let  ||  •  ||  be  an  abbreviation  for 
II  *  l|/€  *  ^  what  follows  x  will  always  denote  a  real  and  z  a  complex  variable}  upper  case 
letters  will  be  used  for  functions  of  x  and  lower  case  for  functions  of  z . 

We  begin  with  a  real  function  F(x)  that  is  continuous  on  It  and  with  a  pair  of  fixed 
integers  m,n  >  0.  Let  Vmn  be  the  set  of  rational  functions  of  type  (m,  n)  with  real 
coefficients,  and  let  /?*(x)  denote  the  best  approximation  to  F  on  /*  out  of  Vmn  .  (  R * 
exists  and  is  unique;  see  [2],  [6],  or  [31].)  For  any  finite  M  >  0,  F  possesses  a  partial 
Chebyshev  expansion 


/ 

F(x)  =  Fm{x)  +  Gm(x)  =  akTk(x/()  -t-  CM(x),  (1.1) 

k=Q 

with  Tk  denoting  the  *  th  Chebyshev  polynomial,  where  the  prime  indicates  that  the 
term  with  k  —  0  should  be  multiplied  by  g  .  Here  ak  is  defined  by  the  inner  product 
[6,  p.  117] 

ak  =  ll  n*mx/c)~  -  . 

.  *  J'f  Vt2  -  X2 

Our  fundamental  transplantation  is  the  map 

x(z)  =  cRez=  it(z  +  z-1),  (1.2) 


a  bijection  of  <1D+  onto  It ,  which  Tor  x  6  /« ,  z  £  DD  leads  to  the  formula 

Tk{x/i)  =  X-(zk  +  z~k).  (1.3) 

In  particular,  let  us  set  o_fc  =  ak  and  define 


Then 


!m{z)  =  Ukzk> 


/+( *)  =  E  akzk, 

*=m— n+l 


rw  = 


m— n 

yi  if  m  >  n, 

k~n—m 

n  — m—  l 

—  ^  0*2*  if  m  <  n. 

fc=xm— n+l 


Fm(x)  =  =  i[/+(z)  +  /+(*-')  +  /•(,)]. 


(1.4a) 


(1.46) 


(1.5) 
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The  idea  of  the  CF  method  is  to  first  approximate  the  analytic  function  /+  on  dD 
by  considering  an  infinite-dimensional  space  Vmn  in  which  a  best  approximation  can  be 
found  exactly,  then  derive  from  this  a  near-best  approximation  Rc *  to  F  on  Ie .  Let 
Vmn  be  the  set  of  functions  that  can  be  written  with  real  coefficients  in  the  form 


m  j  n 

?(*)=  E 

=  — oo  /  Ac=0 


where  the  terms  of  negative  degree  in  the  numerator  converge  to  a  bounded  analytic 
function  in  \z\  >  1  and  the  denominator  has  no  zeros  in  D\JdD  .  Let  H  be  the  real 
symmetric  Hankel  matrix 

/  am— n+ 1  Gm— n+2  ’  *  '  aM  \ 

I  flm-n+2  ] 


(A  Hankel  matrix  is  a  matrix  with  =  aj*+y|  .)  Let 

H  =  UKUt 

be  a  real  orthogonal  eigenvalue  decomposition  of  II  —  i.e.  U  ,  A  are  square  real  matrices 
with  A  =  diag(Xi, . . . ,  \m >n-m)  and  U  is  orthogonal  (  UT  =  U~x  ).  We  assume  the 
eigenvalues  are  ordered  by  absolute  magnitude:  [X^  >  |X2|  >  •••  >  |X/„f+n_m| .  (If  FI 
had  complex  coefficients  it  would  become  clear  that  a  singular  value  dccompositon  rather 
than  an  eigenvalue  decomposition  is  most  appropriate,  but  since  //  is  real  symmetric,  the 
two  arc  the  same  here  except  Tor  the  possibility  of  negative  signs.)  Let  X  abbreviate  Xn+i 
and  let  (tti , . . um +n_m)r  be  the  corresponding  right  eigenvector,  namely  the  (n  +  I)st 
column  of  U  .  The  following  result  was  proved  by  Caratheodory  and  Fcjer  in  1911  for  the 
polynomial  case  n  =  0  [5],  extended  to  rational  approximation  by  Takagi  in  1924  [37], 
and  generalized  further  by  Adamjan,  Arov,  and  Krcin  in  1971  (1).  A  presentation  and 
partial  proof  can  be  found  in  [4lj.  A  full  discussion  of  degenerate  cases  will  be  given  in 
[19|. 

Theorem  1.  /+  has  a  unique  best  approximation  r *  on  dD  out  of  Vmn  >  which 
is  given  by 

f+  -  r%  =  b,  (1.6a) 


where 


The  error  is 


b(z)  =  \zr 


•«i  +•••  +uM+n-mzM+n-m-1 


UM+n-m  +  •••  +  «!  * 


ll/+-?*||aD  =  |X|, 


M  +  n  —  m—  1  * 


and  the  error  curve  (/+  —  r*)(dD)  is  a  perfect  circle  about  the  origin  whose  winding 
number  is  m  +  n  +  l  if  |Xn|  >  |X|  >  |Xn+2| .  | 


t 


The  function  b  is  X  times  a  finite  Blaschke  product,  which  is  why  /+  —  f*  maps 
3D  onto  a  circle,  and  tho  optimality  of  r*  for  the  complex  approximation  problem  can 
be  seen  to  follow  from  this  by  Rouche’s  theorem.  Now  this  optimality  is  not  of  use  to  us. 
However,  let  us  transplant  to  l€  by  defining 


M*)  =  (*)+?*(*  l)  +  /°(2)]. 


(1.7a) 


or  by  (1.5)  and  (1.6a) , 

-  b(z)  ~  M*"1)]- 

Then  by  (1.5)  again, 

Fm(x)  -  R(x)  =  Rc  b(z), 


(1.7  b) 
(1.8) 


and  if  b  has  winding  number  m  +  n  -f  1  ,  it  follows  that  Fm  —  R  equioscillates  on  ft  at 
m  +  n  +  2  points  £  =  x0  >  x\  >  •  •  •  >  xm+n+t  =  — £  : 


\\F*-R\U.  =N, 

(Fm~R)(x3)  =  (— 1)JX. 


(1.9) 


If  R  belonged  to  Vmn  ,  this  equioscillation  would  imply  R  =  R *  and  |\|  =  E*  ,  and  we 
would  have  solved  our  original  approximation  problem  (for  Fm  )  exactly.  Unfortunately, 
this  is  in  general  not  the  case.  (The  main  exception  occurs  when  M  =  m  +  1  ,  and 
this  gives  rise  to  some  of  the  examples  of  Achieser,  Talbot,  and  others  mentioned  in  the 
Introduction.)  But  the  key  to  the  CF  method  is  that  for  smooth  functions  F  ,  R  turns 
out  to  be  very  close  to  Kmn  . 

Let  q  =  *kzk  denote  the  normalized  denominator  of  r*  —  the  polynomial  of 

degree  dq  <  n  with  constant  term  =  1  whose  zeros  are  the  finite  poles  of  r*  lying 
outside  3D  .  Define 

Q(z)  =  <i(zMz~1)/t,  (1.10) 

where  r  is  the  scalar  q(i)q(—i) ,  inserted  to  make  Q(x)  have  constant  term  1  .  Note  that 
Q[x)  =  \q{z)\2jr  >  0  on  /e .  Now  since  ?*  €  Vmn  , 

f+(z)  —  b(z)  —  r*(z)  =  Q(zjn~dq)  as  z  — ►  oo, 


hence  since  dq  <  n  and  f m  —  /+  =  0(zrn^n) , 

Jm{z)  —  b(z)  —  0(zrn~dq)  as  z  — ►  oo. 

Let  us  consider  the  Laurent  series  with  respect  to  3D  of  the  product 
R(x)Q(x)  =  ^[/m(z)  -  b(z)  -  b(z~l)]q(z)q(z~l)/T. 


(1.11) 


(1.12) 
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By  (1.11)  and  the  definition  of  q ,  [/m(z)  —  b(z)]q(z)q(z~l)/r  must  be  analytic  outside 
dD  except  for  a  pole  of  order  at  most  m  at  oo  ,  and  therefore  all  terms  of  order  greater 
than  m  in  the  Laurent  series  of  (1.12)  are  due  to  b[z~x  )q(z)q(z~l)/r  .  By  symmetry,  all 
terms  of  order  less  than  — m  are  due  to  b(z)q(z)q(z~{)/r .  Hence,  if  we  define 


M 


b(z)  =  b(z)q{z)q(z  *)/r  =  ^  &***> 


k=—oo 


b  =  ^2  bkzk>  b  (z)  =  b(z)  ~  b  (2)> 

Jk=s— m 

then  the  function 

p‘(x) =  “**(*'“  **(*”*) 
is  a  polynomial  of  degree  m  in  x .  If  we  further  set 


"  9$  ’  bR(i1  =  5 


fr^(z)  -f*  6^(2:  *) 


then  we  obtain 


hence 


/?(*)£(*)  =  «t7Q(*)  +  0(Tm+,(*)). 


(1.13a) 

(1.136) 

(1.14a) 

(1.146) 

(1.15) 


We  will  call  R\*  G  the  type  l  or  Macbly  type  CF  approximation  of  F  ,  because 
as  in  Machly’s  generalization  of  Fade  approximation  (cf.  [6],  p.  1 18  and  [17]),  truncation  of 
higher-order  terms  in  R  is  done  after  multiplying  through  by  the  denominator  Q  .  There 
is  a  second,  probably  superior  way  to  truncate  R ,  namely  by  using  a  Chcbyshev-Pade 
kind  of  approximation  with  fixed  denominator  Q  .  That  is,  one  may  take  the  type  2  or 
Fade  type  CF  approximant  as 

Ri'M  =  9$. 

with  P2  defined  by  the  condition 


R(x)  =  Rf/(x)  +  0(Tm+i(x)). 


(1.166) 


One  could  go  further,  in  complete  analogy  with  the  Chcbyshev-Pade  approximation  defined 
by  Cragg,  and  define  a  third  type  of  CF  approximation  by  permitting  the  denominator  of 
Rc?  to  be  free  as  well  as  the  numerator.  However,  one  might  then  end  up  with  a  fraction 
having  a  pole  on  It .  For  this  reason,  and  on  the  basis  of  numerical  experiments  and  the 
analogy  with  the  Chcbyshev-Pade  situation,  we  believe  that  R*J  is  the  best  of  these  three 
possibilities,  and  from  now  on  we  will  drop  the  subscripts  and  assume  Rc *  =  R\*  . 
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To  obtain  the  polynomial  P  (  =  P2  )  satisfying  (1.16),  one  proceeds  as  follows.  Let 
Ck  ,  0k,  and  7*  denote  the  A:  th  Chebyshev  coefficients  of  R(x) ,  P(x) ,  and  1  /Q(x), 
respectively.  Then  the  coefficients  {0k}  satisfy  the  Toeplitz  system  of  equations 

(  10 

li 

V  72m 

Since  \/Q(x)  >  0  on  I€  ,  the  infinite  symmetric  Toeplitz  matrix  (7|t-y|)^5=-00  is  known 
to  be  positive  definite  [18],  and  hence  the  principal  submatrix  appearing  in  (1.17)  is  positive 
definite  also,  hence  nonsingular.  Moreover,  since  this  submatrix  is  symmetric  about  its 
anti-diagonal,  a  symmetric  right- hand-side  (ik)  leads  to  a  symmetric  solution  (/?*) ,  as 
indicated  in  (1.17).  Indeed,  in  practice  one  may  reduce  (1.17)  to  a  system  of  size  m  +  1 
instead  of  2m  +  l  .  Consequently,  P  is  always  well  defined  by  (1.16). 

In  summary,  here  is  the  real  rational  (type  2)  CF  approximation  method.  We  have 
indicated  lour  points  at  which  a  numerical  implementation  can  naturally  be  based  upon  the 
Fast  Fourier  Transform.  For  further  information  on  uses  of  the  FFT  in  complex  analysis, 
see  [22].  The  FFT  method  indicated  for  the  polynomial  factorization  of  Step  4  is  that 
proposed  in  §3.2  of  [22] ;  see  also  [41]. 

Step  1.  Given  F ,  find  its  Chebyshev  coefficients  do,  •••»%  for  some  large  M  (FFT), 

Step  2.  Construct  the  Ifnnkel  matrix  //  and  find  its  n  +  1  st  eigenvalue  (in  absolute 
value)  and  eigenvector . 

Step  3.  Find  the  Laurent  series  on  the  circle  of  the  Blaschkc  product  b(z)  defined 
by  (l.6b)  (FFT).  Subtract  this  plus  its  conjugate  from  /m{z)  to  obtain  the 
Chebyshev  coefficients  {c*}  for  R(x)  by  (l.7b). 

Step  4.  Factor  the  denominator  of  (l.Gb)  to  obtain  the  polynomial  q(z)  and  construct 
Q(x)  from  (1. 10)  (FFT).  Find  the  Chebyshev  coefficients  {7^}  for  i/Q(x)  (FFT), 

Step  5.  Determine  the  polynomial  P{x)  satisfying  (L16)  by  solving  (1.17),  and  define 
R'f  =  P/Q . 

Step  6.  To  get  a  bound  on  Ecf  —  E*  ,  examine  how  close  the  error  curve  of  Rr *  comes 
to  cquioscillating . 


Remark.  This  somewhat  obscure  construction  of  Rc *  can  be  made  much  more 
transparent  in  the  case  m  >  n  .  The  theory  of  complex  CF  approximation  shows  that  f* 
in  Theorem  l  is  close  to  Vmn  ,  and  for  m  >  n,  (l.7a)  then  implies  that  R  is  close  to 
Vmn  also. 


2.  ASYMPTOTIC  RESULTS  FOR  SMALL  INTERVALS 


The  basis  of  our  results  Tor  small  intervals  is  the  theory  worked  out  in  [41]  for  complex 
approximation  on  small  disks.  For  these  results  a  normality  assumption  is  needed.  Let  F 
be  given,  and  let  M  >  3m  +  2n  4*  2  be  a  fixed  integer. 


Assumption  A.  The  M  th  derivative  of  F(x)  exists  and  is  Lipschitz  continuous  at 
x  =  0  .  Moreover ,  if  F(x)  —  +  0(xM+l) ,  with  a*  =  0  for  k  <  0 ,  then 


\  '  1  '  — 2  / 


Equivalently,  the  Fade  approximation  of  F  of  type  (m,n)  has  a  full  n  finite  poles — see 
[23],  Thms.  7.5c  f  and  [16],  §3. 


The  nonvanishing  of  this  Hankel  determinant  is  a  standard  assumption  that  appears 
also,  for  example,  in  [12],  [43],  and  [31,  p.  170].  For  many  functions,  including  ex  ,  it  is 
satisfied  for  all  (m,  n) . 

Here  is  the  main  result  from  [41]  that  we  require: 

Theorem  2.  Let  F  satisfy  Assumption  A,  and  for  each  e  >  0 ,  let  r*  be  the 
extended  best  approximation  in  Vrnrx  of  the  function  /+  defined  by  (1,4a)  .  Then  for 
ail  sufficiently  small  t  ,  b  =  /+  —  r*  has  winding  number  exactly  m  +  n  +  1  on  dD ,  it 
approximates  a  monomial  according  to 

(/+  -  r'){z)  =  \;m+n  +  ,(l  +  0{c))  (2.1) 


uniformly  on  OD  ,  and  its  Laurent  coefficients  on  OD  satisfy 

bk  =  [0(e)]2m+2n+2-fc  V  k  <  Til  +  71  +  1  (2.2a) 


uniformly  in  k  .  In  addition ,  the  coefficients  of  the  denominator  q  of  ?*  satisfy 

rk  —  0(ek)  (0  <  k  <  n).  (2.26) 

Proof.  It  can  be  seen  that  as  e  — ♦  0  ,  {a*}  and  {ak}  are  related  by 

ak  =  a|fc|2'-|fc|t|fcl  +0(f|fc|+1)  -  M  <  k  <  M 

(cf.  (21),  Lemma  3.3),  hence  since  afe  =  0  for  k  <  0, 

ak  =  ak2l-kek  +  0(ffc+l)  -  M  <  k  <  M.  (2.3) 

Now  if  the  0(cfc+l)  term  were  zero,  the  extended  approximation  problem  for  /+  would 
bo  that  of  approximating  ]C*Lo  a*2,_fc(r£)fc  ,  and  for  this  the  theory  of  [41]  applies.  From 
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Lemmas  1.1,  1.3  and  4.4  of  [41]  one  would  obtain  (2.26),  (2.2u) ,  and  (2.1),  respectively, 
on  the  basis  of  Assumption  A.  In  fact  the  term  0(cfc+1)  is  not  zero,  but  it  is  of  size  G[t) 
relative  to  the  terms  just  considered.  This  is  enough  to  make  the  arguments  of  [41]  still 
go  through;  we  omit  the  details.  In  particular,  Assumption  A  and  (2.3)  imply  that  for  all 
sufficiently  small  e ,  the  corresponding  Hankcl  matrix  made  up  of  coefficients  a ^  will  also 
have  nonzero  determinant.  | 

We  will  need  a  lemma  on  the  behavior  of  the  denominators  of  Rc *  and  /?*  .  Let  us 
write 

Rcf  =  P/Q,  «*  =  P'/Q\ 

where  P  and  P *  are  polynomials  of  degree  at  most  m  ,  and  Q  and  Q *  arc  polynomials 


of  degree  at  most  n  with  constant  terms  1  . 

Lemma  3  (cf.  [ 41 ],  Lemma  6.1).  As  e  — ►  0  , 

<3=1+  O(e), 

(2,1a) 

<3*  =  i  +  o(c), 

(2.46) 

and  hence 

QQ’  =  l  +  O(c) 

(2.4c) 

uniformly  on  fc  . 

Proof.  Equation  (2.1a)  follows  from  (2.26).  Equation  (2.16)  is  a  corollary  of  the 
known  fact  [7,13]  that  as  f  — *  0  ,  R*  approaches  coefficientwise  the  Pade  approximant 
Rp  G  Vmn  i  whose  normalized  denominator,  having  coefficients  independent  of  r,  ob¬ 
viously  satisfies  Qp  —  1  4*  O(c)  uniformly  on  /«  as  c  — ►  0  .  (In  fact  [7]  and  [43]  show 
R *  — ►  Rp  for  [0,t],  not  [ — 6 ,  ^ ]  .  However,  C.  Chui  of  [7]  assures  us  that  the  result  also 
holds  for  the  latter  problem  (private  communication,  September  1981).)  | 

We  now  show  that  from  Theorem  2  it  follows  that  F  —  Rcf  equiosci Hates  at  m4-n  +  2 
points  on  It  to  within  0(t3m+2n+3) .  Let  {zj}  be  the  set  of  (2m  +  2n  4-  2)  th  roots  of 
unity  on  dD+  ,  and  let  {xj}  be  the  corresponding  set  of  Chebyshev  abscissae  for  /«  , 
xj  =  e  R ezj  . 

Theorem  4.  Let  F  satisfy  Asssmption  A ,  and  for  each  c  >  0,  Jet  Rc *  be  its 
CF  approximation  out  of  Vmn  .  Then  for  all  sufficiently  small  c ,  there  is  a  set  of  points 
C  =  x0  >  X\  >  •••  >  xm+n+t  sss  -I  satisfying 

xj  =  xj(l  +  0(c))  (2.5) 

at  which 

||F  -  Rcf  ||  -  sgn  X(-l)>(F  -  Rcf)(xj)  =  0(e3m+2n+3).  (2.6) 

Proof.  Lot  {z^  be  the  sot  of  m  +  n  +  2  points  on  dD+  at  which  /+  —  f*  is 
real,  numbered  in  Counterclockwise  order  from  zq  =  1  to  zm+n+ i  =  —1 ,  and  take 
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Xj  =  e  Rezj  .  The  bound  (2.5)  Follows  from  (2.1).  By  (1.9),  Fm  —  it  exactly  equioscillates 
on  the  set  {x;}  with  error  |X| ,  so  to  prove  (2.6),  it  is  enough  to  show 


||F  -  Fm ||  =  o(£3m+2n+3) 

(2.7) 

l|/e-/?^H  =  c»(£3m+2"+3). 

(2.8) 

The  first  bound  Follows  From  the  Lipschitz  continuity  statement  of  Assumption  A.  To 
establish  (2.8),  wc  observe  that  from  (1.13a)  and  (2.2a,  6) ,  we  have 

bk  =  [0(£)]2m+2n+2-fc  V  k  <  m  +  1. 

This  implies  by  (1.136) , 

bR  =  0(e3m+2n+3). 

From  this  and  (1.14)  and  (2.4a)  it  follows  that  the  analog  of  (2.8)  holds  for  type  1  CF 
approximation: 

||ff  -  R\f  ||  =  <9(£3m+2n+3).  (2.9) 

*  ft 

Now  6  has  terms  only  of  degree  <  — m  —  1  ,  each  of  order  0(e3m+2n+3) ,  and  from 
(1.10)  and  (2.26)  it  is  straightforward  to  see  further  that  the  Chcbyshev  coefficients  of 
l /Q{x)  satisfy 

b*l  =  Io(c))|fc|  v  k. 

From  these  facts  and  (1.14)  it  follows  that  the  degree- m  part  of  the  Chebyshcv  series  of 
it  —  !t\ *  has  magnitude  0(£3m+2n'f  4) .  In  other  words,  if  the  j3  -coefficient  vector  For  Pi 
is  inserted  in  the  left  hand  side  of  (1.17),  then  that  system  is  satisfied  up  to  an  error  of 
magnitude  0(c3m-f2n“M) .  But  as  c  — ►  0  ,  the  matrix  in  (1.17)  approaches  the  identity, 
so  it  Follows  that  the  coefficients  of  Pj  agree  with  those  of  P  up  to  (7(t3m^2n+1) .  By  a 
final  application  of  (2.1a) ,  this  implies 

|| R\f  -  Re} ||  =  o(£3m+2n+4) 

and  with  (2.9)  this  yields  (2.8).  | 

Thus  the  error  curve  of  Rc ?  equioscillates  up  to  0(€3m+2n+3) .  By  the  de  la  Vall6c 
Poussin  theorem  for  rational  approximation  [31,  Thm.  98j,  this  implies 

Corollary.  As  c  — ►  0 

||P  -  Pc/||  -  ||P  -  R'\\  =  tf(c3m+2n+3) 


and 

|X|  -  ||F  -  P*||  =  0(c3m+2n+3).  | 
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Wo  wish  to  show  further  that  \\RC*  —  /?*||  =  0(€3m+2n+3) .  Now  by  definition 
\\F  —  R* ||  <  \\F  -  /ic^|| ,  so  (2.5)  implies  that  at  the  m  +  n  +  2  points  {x;}  ,  Ref  —  R * 
must  satisfy  an  alternating  sequence  of  constraints 

-sgnX(— 1  y(lic'  -  /?*)( Xj)  <  v,  (2.10) 

where  7/  =  0(£3m+2n+3)  Does  this  imply  that  Rcf  —  Rm  is  small?  In  fact  it  does,  and  this 
question  was  taken  up  directly  in  the  paper  on  rational  approximation  on  small  intervals 
by  Machly  and  Witzgall  [30,  Lemma  4.6].  however,  for  our  needs  (2.10)  can  be  reduced 
to  a  similar  set  of  constraints  on  a  polynomial  instead  of  a  rational  function,  which  will 
be  easier  to  deal  with.  Let  us  write 

Rd  _  t  _  H  -  ggl  ~ 

Q  Q *  QQ"  ’ 

and  let  S  denote  PQ*  —  P*Q  ,  a  polynomial  of  degree  at  most  m  +  n  .  Then  by  (2.4c) , 

Rc*  -  R‘  =  S{1  +  0{c))  (2.11) 

uniformly  on  fe  as  e  — ►  0  .  Therefore  (2.10)  leads  to  the  sequence  of  m  +  n  +  2  constraints 

-sgn\(-l),S(iJ)  <  rf  (2.12) 

for  some  new  r/  =  0(f3m+2n+3) .  We  want  to  deduce  that  ||S||  =  0(<3m+2n+3) . 

This  is  a  commonly  occurring  problem  in  approximation  theoretic  proofs.  If  {xj} 
were  a  fixed  set  of  points  (i.c.  not  dependent  on  e ),  then  the  argument  that  would  be 
required  is  the  key  step  in  proofs  of  strong  uniqueness  or  Lipschitz  continuity  for  poly¬ 
nomial  Chebyshov  approximation.  Essentially  the  same  reasoning  for  this  has  appeared 
in  (at  least)  papers  of  Freud  [15],  Machly  and  Witzgall  [29],  and  Cline  [8];  Machly  and 
Witzgall  even  give  a  figure  illustrating  (2.12)  graphically.  For  a  general  discussion  see  [6]. 
In  our  application  the  near-alternation  points  are  not  fixed,  but  by  (2.5)  they  arc  close 
to  Chcbyshev  abscissae  for  small  c  ,  hence  uniformly  separated  from  each  other.  This 
uniform  separation  is  what  is  needed  to  make  the  argument  go  through,  and  the  same  is 
the  case  for  applications  to  strong  uniqueness  and  Lipschitz  continuity  [11,24], 

Lemma  5  (cf.  [*2lj,  Thm  2.1).  Let  S(x)  be  a  real  polynomial  of  degree  at  most  p 
on  I€  .  Suppose  there  exist  p  +  2  points  e  >  xq  >  x\  >  •••  >  xM+i  >  —  €  at  which 


(-l)»S(*y)  <  V 

(2.13) 

for  some  rj  >  0  ,  and  suppose  that  Xj  =  c  cos  <\>j  with 

6  3*  ^  6 

J  H  +  1  “  p(2n  +  1) 

(2-14) 

for  some  6  <  1  .  Then 

IISII  < 

(2.15) 
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Proof.  S(t  cos0)  is  a  trigonometric  polynomial  of  degree  p  in  <j> ,  so  by  Bernstein's 
inequality  [6,  p.  91],  one  has 

<  rfISII. 

If  {zj}  are  the  Chebyshcv  points  xj  =  ecos(jir/(fi+  1)) ,  then  this  bound  together  with 
(2.13)  and  (2.14)  implies 

(-l)'S(xf)  <  k  (2.16) 


where 


K  =  Tf  + 


'  2jt  +  l 

Now  according  to  a  computation  of  Cline  in  [8],  §4  ,  (2.16)  implies 

11*511  <  (2^+l)«. 


(2.17) 


Therefore 

l|S||  <  (2/i  +  l){l  +  =  (2/x  +  1)„  +  6\\S\\, 

hence  (2.15).  | 

Lemma  5  provides  all  that  is  needed  to  prove  our  first  main  theorem: 

Theorem  6.  Let  F  satisfy  Assumption  A,  and  for  each  e  >  0  let  Rc?  and  R*  be 
its  CF  and  Chebyshcv  approximations  in  Vmn  .  Then  as  e  — *•  0  , 

|| Rcf  -  /T||  =  0(c3m+2n+3). 

Proof.  Applying  Lemma  5  to  (2.12)  with  p  =  m  +  n  gives 

||5||  =  0(c3m+2n+3). 

The  result  follows  from  (2.11).  | 

Together  with  (2.8)  Theorem  6  implies 

||JT-  R\\  *  0(c3m+2n+3).  (2.18) 

Wc  can  interpret  this  as  a  statement  about  the  geometry  of  optimal  error  curves,  analogous 
to  the  theorems  in  [-10]  and  [41]  showing  that  error  curves  in  complex  Chebyshcv  ap¬ 
proximation  arc  close  to  perfect  circles: 

Theorem  7.  Let  F  satisfy  Assumption  A.  Then  for  all  sufficiently  small  c ,  there 
exists  a  function  b(z)  that  is  a  constant  multiple  of  a  Unite  Dlaschkc  product ,  is  analytic 
in  1  <  |z|  <  oo  except  for  at  most  n  poles ,  and  has  winding  number  m-f  n+  1  on  dD , 
satisfying 

||(F  -  /f*)(z)  -  Rc  6(z)||  =  0(c3m+2n+3). 
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Proof.  Follows  from  (1.8),  (2.7),  and  (2.18).  | 

Theorems  6  and  7  are  the  extensions  to  rational  approximation  of  results  given  in 
Theorem  3.4  of  [21].  Theorem  3.5  of  that  paper  also  proves  analogous  estimates  for  m  — ► 
oo  on  a  fixed  interval,  but  we  have  not  extended  these  results. 

If  F  satisifes  Assumption  A,  then  ||F  —  R*\\  has  size  0((m+n+i)  but  not  o(tm+n+l) 
as  £  — ►  0  .  Relative  to  this  scale,  therefore,  our  results  have  strength  0(e2m+n+2).  It 
appears  that  these  orders  are  best  possible,  except  that  in  the  case  n  =  0 ,  a  certain 
“bonus”  cancellation  makes  it  possible  to  increase  3m  +  3  to  3m  +  4  in  Theorems  6  and 
7  (but  not  in  the  estimate  on  |X|  in  the  corollary  to  Theorem  4).  See  the  proof  of  Lemma 
3.2  in  [21]  for  details. 
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3.  NUMERICAL  EXAMPLES 

The  CF  method  is  not  difficult  to  implement  numerically,  and  the  techniques  we  have 
used  are  described  in  §7  of  [41].  In  outline,  we  rely  heavily  on  the  Fast  Fourier  Transform 
as  indicated  here  at  the  end  of  §1  ,  and  use  EISPACK  routines  based  on  Sturm  sequencing 
for  the  eigenanalysis  of  Step  2.  The  bottleneck  is  the  eigenvalue  computation,  which  takes 
time  0(M3) ,  for  unfortunately  no  way  is  known  to  take  advantage  of  the  Ilankel  structure 
of  H  .  Ideally  one  wants  M  to  be  large  enough  so  that  the  Chebyshev  coefficients  a *  for 
k  >  M  are  negligible,  hence  Fm  «  F ,  For  most  of  the  examples  considered  below,  such 
as  those  involving  ex  ,  this  is  achieved  with  Af  =  35  ,  leading  to  computation  times  on  the 
order  of  0.1  secs,  on  our  IBM  370/168.  With  F(x)  =  |z| ,  on  the  other  hand,  M  =  120 
is  only  barely  large  enough  to  get  approximations  with  m,  n  <  2  accurately,  and  the 
computation  time  increases  to  2  secs.  Thus  the  CF  method  is  not  only  more  accurate,  but 
also  much  faster  if  the  function  to  be  approximated  is  smooth.  For  certain  high- precision 
numbers  below  we  have  resorted  to  quadruple  precision. 

In  general  the  CF  method  will  yield  an  approximation  satisfying 

|X|  «  E\ 


and  by  the  de  la  Vallee  Poussin  result  one  has 


<  nef, 


(3.1) 


where 

i„  =  "in  | (F  -  Rcf){xj)\, 

J 


with  the  minimum  taken  over  that  set  of  m  +  n  +  2  nearly  alternating  points  which 
maximizes  its  value.  In  our  experiments  Rc *  and  X  were  computed  to  close  to  machine 
precision,  and  by  means  of  a  minimization  routine  (FMIN,  by  Richard  Brent),  Ec *  and 
^min  wcre  also  found  to  this  accuracy  (Step  6).  The  quantity  not  so  precisely  known 
is  E *  ,  for  we  do  not  have  a  high-accuracy  rational  Chebyshev  approximation  routine  at 
hand.  Therefore  in  what  follows  we  report  |\|  rather  than  E*  . 

As  a  first  example,  let  us  give  more  details  related  to  ex  .  Table  2  shows  the  eigenvalue 
|X|  in  approximation  on  [—1,  l]  for  0  <  m,  n  <  3.  Each  digit  known  to  agree  with  the 
corresponding  digit  of  E*  after  both  are  rounded  has  been  underlined.  In  most  cases 
this  knowledge  is  based  on  (3.1).  The  agreement  is  excellent,  and  we  believe  that  it  will 
get  steadily  better  as  m,  n  — oo  in  any  fashion.  (The  table  leaves  some  doubt  as  to 
whether  this  is  true  for,  say,  m  =  0  and  n  — *  oo  ,  but  further  experiments  show  that 
by  type  (0,9),  five  significant  digits  can  be  underlined,  and  the  number  is  growing  with 
n .)  This  example  relates  to  a  conjecture  of  Mcinardus  [31,  p.  168]  which  proposes  that 
as  m,n  — *■  oo  , 


/5* 

y  mn 


(«*)« 


m!  n! 


2m+n(rn  +  n)!  (m  +  n  -f  1)! 


(l+o(l))- 
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(3.2) 


m  *  0 

m  =  1 

m  =  2 

TO  =  3 

=  0 

1.1961(-0) 

2. 787994(-l) 

4.5017  38776 (-2) 

5.5283 

70108  71194(5-3) 

=  1 

2. 1724  C-l) 

2.096982(-2) 

1.7890  66755<-3) 

1.3461 

23369  20018(-4) 

=  2 

3 ■ 5288 (—2 ) 

1.677017(-3) 

8.6899  91075 (-5) 

4.3991 

63371  96896 (-6) 

=  3 

_4.  5235(-3) 

1.239861(-4) 

4.2766  46704(-6) 

1.5506 

69053  97117 (-7) 

Table  2.  The  eigenvalue  \  A  |  in  CF  approximation  to  eX  on  (jl,l  j 
for  various  (m,n)  .  Underlined  digits  are  known  to  agree  with 
corresponding  digits  of  E*  ,  after  rounding.  Doubly  underlined 
digits  agree  with  the  conjectured  limit  formula  (3.2)  of  Meinardus. 


At  this  writing  the  conjecture  is  unproved,  but  it  is  known  to  be  valid  up  to  a  constant 
'  factor  [4,33]  of  less  than  40.  The  double  undcrlinings  in  Tabic  2  mark  digits  where  (3.2) 
agrees  with  |X| ,  and  it  is  evident  that  |X|  is  much  closer  to  E *  than  (3.2)  is.  This 
suggests  that  it  might  be  possible  to  resolve  the  Meinardus  conjecture  by  means  of  the  CF 
method.  (Chebyshcv-Pade  approximation,  it  turns  out,  is  not  strong  enough  (17].)  Two 
key  proofs  would  be  enough  to  settle  the  issue: 

0)  IM  =  £mn(l  +  °(1))  as  w»,  n  —  OO, 

(ii)  |X|  satisfies  (3.2)  as  m,  n  — *  oo  . 

Unfortunately,  these  claims  are  not  at  all  easy  to  establish,  and  the  Meinardus  conjecture 
will  probably  be  proved  true  before  long  by  some  simpler  technique.  Nevertheless,  it  would 
be  very  interesting  to  know  that  asymptotically,  best  approximation  errors  agree  with  the 
eigenvalues  of  an  infinite  Hankcl  matrix  of  Chebyshev  coefficients. 

Incidentally,  we  conjecture  that  for  the  problem  of  complex  approximation  of  e*  on 
the  unit  disk,  formula  (3.2)  holds  with  the  factor  2m+n  removed.  The  best  approximation 
errors  computed  in  [41],  Tabic  3  by  the  complex  CF  method  show  that  the  agreement  of 
this  conjecture  with  exact  best  approximation  errors  for  m,  n  <  3  is  about  as  close  as  in 
the  real  case.  The  disappearance  of  the  power  of  2  is  natural  in  the  light  of  (2.3). 
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Tabic  3  shows  some  results  of  the  CF  method  for  functions  besides  ex  .  In  each 
case  |Xj  is  given  for  approximation  of  type  (0,1),  (1,1)  and  (2,1)  on  j— 1,1),  with 
underlinirigs  as  before.  Evidently  CF  approximation  really  works,  even  for  as  low  a  degree 
as  (l,  l) .  It  is  apparent  that  it  performs  relatively  poorly  for  |x| ,  and  this  is  not  surprising 
in  view  of  the  nondirfcrcntiabilily  of  this  function. 

Now  let  us  confirm  that  the  asymptotic  orders  of  accuracy  predicted  in  Section  2 
arc  valid,  and  sharp.  Table  4  shows  (X j  and  FZcf  —  £^in  0*0  approximation  of  e% 
on  I€  for  €  =  With  each  factor  of  2  reduction  in  6,  we  expect  |X(  to 

decrease  by  approximately  2m'fn+1  =  8  ,  and  Ec *  ~  Ec^in  to  decrease  by  approximately 


23m-4-2n-*-3  =_=  256 

The  table  confirms 

these  predictions.  Such  asymptotic  behavior, 

F(x) 

Ul  :  (0,1) 

Ul:  (1,1) 

r-l 

CnI 

r< 

Ui 

4. 483(-l) 

4. 4827(-l) 

1.1359(-1) 

6 

X 

5. 397 (-1) 

5.3970(-l) 

1.9257(-1) 

/l.  1  x 

2.238(-l) 

1.633K-2) 

2.9709 (-3) 

arctan(x) 

8. 312 (-1) 

4. 7889 (-2) 

4. 7889(-2) 

i/r(x+i) 

4.041(-1) 

1. 1955(-1) 

2.1045  75498(-2) 

X 

e 

2.172(-1) 

2.0970(-2) 

1.7890  66755(-3) 

i  x+3 

log  £ 

1. 598(-l) 

8.6079  41336 (-4) 

4.9591  156l392(-5) 

Table  3.  The  eigenvalue  |Xj  in  CF  approximation  of  types 
(0,1),  (1,1)  and  (2,1)  to  various  functions  on  01,1^  • 

Underlined  digits  are  known  to  agree  with  corresponding  digits 
of  E*  , 
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e 

N 

Ratio 

Ecf  -  Ecf 
mm 

Ratio 

1 

2 . 097(-2) 

2.03  (— 6 ) 

1/2 

2 . 605 (—3) 

8.04 

9 . 18 (—9 ) 

221 

1/4 

3.255 (-4) 

8.00 

3. 7  3 (—11) 

246 

1/8 

4.069(-5) 

8.00 

1.47 (-13) 

253 

1/16 

5.086C-6) 

8.00 

5. 77 (-16) 

255 

Table  4. 

Type  (1,1) 

approximation 

of  e  on 

T-e.eJ  for 

various 

e  . 

however,  is  dependent  on  the  smoothness  of  F  .  With  F(x)  =  |x| ,  for  example,  cutting 
c  in  half  is  equivalent  to  dividing  F  by  2  ,  and  therefore  the  effect  will  be  to  cut  X  and 
Ecf  —  Ec}n\n  exactly  in  half,  no  matter  what  m  and  n  are.  Thus  for  F{x)  =  |x|  there  is 
nothing  to  be  gained  by  shrinking  the  interval,  although  the  CF  method  still  improves  if 
m  and  n  are  increased. 

Our  final  example  is  associated  with  some  further  conjectures  about  asymptotic 
degree  of  approximation.  In  a  paper  of  Cody,  Meinardus,  and  Varga  [9],  the  problem  of 
approximating  e~e  on  the  semi-infinite  interval  [0,  oo)  was  studied.  They  proved  that 
in  rational  approximation  of  type  (0,n)  or  (n,  n) ,  the  error  decreases  geometrically  as 
n  — ►  oo  ,  but  did  not  determine  an  asymptotic  rate  of  decrease.  For  approximation  of  type 
(0,n)  ,  they  gave  numerical  results  that  suggested  the  limiting  behavior 


and  this  equality  was  later  proved  valid  by  A.  Schonhagc  [35].  For  approximation  of  type 
(n,n),  their  numerical  results  reported  for  n  <  14  suggest  to  us 
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(3.3) 


lirf»  (/?;„),/n  =  B« 


9.28...' 


and  a  limit  1/9  has  also  been  conjectured.  But  no  result  of  this  kind  has  been  established. 

Despite  appearances,  this  problem  can  be  approached  by  the  CF  method.  The 
function 


maps  [0,  oo]  bijcctively  onto  (—1,1],  inducing  a  one-to-one  correspondence  between  ra¬ 
tional  functions  It(t)  and  R(x)  in  Vnn  .  Under  this  transplantation  an  equioscillating 
curve  on  one  domain  maps  to  an  equioscillating  curve  on  the  other,  and  as  a  consequence 
it  can  be  shown  that  the  given  problem  is  equivalent  to  the  problem  of  approximating 

F(x)  =  e(l“l)/(x_M) 


on  [—1,1],  which  can  be  treated  by  the  CF  method.  (For  justification  sec  the  theory 
of  [3];  this  transplantation  only  works  when  m  =  n .)  Here  F  is  C°°  on  [—1,1),  but 
analyse  only  on  (—1, 1] ,  so  its  Chcbyshev  series  decreases  faster  than  any  polynomial 
but  not  geometrically.  We  took  M  =  200  ,  leading  to  truncated  terms  |a*|  <  10~23  for 
k  >  M  ,  a^d  calculated  eigenvalues  in  quadruple  precision  on  the  IBM  370/168.  This  gave 
us  \n+i  for  0  <  n  <  18  accurate  to  many  places,  and  for  n  >  5,  |X|  agrees  with  the 
vrd\ k-  V*  reported  by  Cody,  Meinardus  and  Varga  to  all  four  places  that  they  give.  We 
have  little  doubt  that  |\|  =  Z?*(l  4-  o( l ))  as  n  — ►  oo  . 

On  the  basis  of  these  numbers  we  conjecture 


«+l 1 1  _ 

■'nn  J  U 


for  some  constants  A  ^  .656 ,  B  ffs  9.28903  .  The  evidence  is  summarized  in  Table  5. 
Remarkably,  in  the  course  of  this  writing  Schonhage  [36)  has  independently  conjectured 
that  (3.3)  approaches  a  limit 


|(2  -  vS)’  ~ 


1 

9.28547* 


In  fact  he  proves  that  lim  inf(Z?*n)1/n  is  at  least  two-thirds  of  this  value.  The  closeness 
but  inequality  of  these  two  conjectured  limits  is  striking. 
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Ratio 


n 

Ratio 

p  =  |x  A  , , ! 

Richardson  extrapolant 

P_  +  t(p„  -  Pn/7) 
n  in  n/Z 

0 

.560172  (-0) 

1 

.668057  (-1) 

8.38508 

2 

.735558  (-2) 

9.08232 

9.31473 

3 

.799452  (-3) 

9.20078 

4 

.865210  (-4) 

9.23998 

9.29253 

5 

.934574  (-5) 

9.25779 

6 

.100845  (-5) 

9.26740 

9.28961 

7 

.108750  (-6) 

9.27316 

8 

.117227  (-7) 

9.27689 

9.28920 

9 

.126329  (-8) 

9.27944 

10 

.136112  (-9) 

9.28127 

9.28910 

11 

.146631  (-10) 

9.28262 

12 

.157946  (-11) 

9.28364 

9.28905 

13 

.170119  (-12) 

9/28444 

14 

.183217  (-13) 

9.28507 

9.28904 

15 

.197314  (-14) 

9.28558 

16 

.212485  (-15) 

9 . 28600 

9.28904 

17 

.228815  (-16) 

9.28635 

18 

.246392  (-17) 

9.28664 

9.28903 

Table  5,  Eigenvalues  |x|  in  (n,n)  approximation  of  exp(^— ) 
on  for  various  n  .  This  is  equivalent  to  approximation 

of  e  t  on  [O,00)  .  Digits  agreeing  with  the  values  of  E*  given 
by  Cody,  Meinardus ,  and  Varga  [9  ~]  have  been  underlined,  (The 
results  of  \_9  ]  were  given  to  four  places  and  for  n<  14  only.) 

The  ratio  appears  to  approach  a  limiting  value  9.28903.... 
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